Sub Shot-Noise interferometric phase sensitivity 
with Beryllium ions Schrodinger Cat States 
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Interferometry with NOON quantum states can provide unbiased phase estimation with a sensi- 
tivity scaling as A9 ~ 1/Nt given a prior knowledge that the true phase shift 6 lies in the interval 
—n<9<n. The protocol requires a total of Nt — 2 P — 1 particles (unequally) distributed 
among p independent measurements and overcomes basic difficulties present in previously proposed 
approaches. We demonstrate the possibility to obtain a phase sensitivity beating the classical shot- 
noise limit using published probabilities retrieved experimentally for the creation of Schrodinger cat 
quantum states containing up to N = 6 beryllium ions. 



Introduction. Interferometry plays a central role in 
the development of basic science and new technologies. 
Its main goal is to estimate phase shifts generated by the 
interaction of the interferometer with some external per- 
turbation in domains spanning from micro-scales, as in 
the measurement of Casimir forces, to cosmic-scales, as 
in the detection of gravitational waves. The precise esti- 
mation of phases is limited by two quite different sources 
of noise. Classical noise can be created by micro-seismic 
geological activities, temperature fluctuations, poor de- 
tection efficiencies which, in principle, can be arbitrarily 
reduced. A second source of uncertainty is provided by 
the laws of Quantum Mechanics, and cannot be reduced 
beyond the limits imposed by Heisenberg uncertainty re- 
lations and the Cramer-Rao lower bound [l[ . Interferom- 
etry with uncorrelated particles allows phase estimations 
with sensitivity bounded by the standard quantum limit 
(shot-noise) 
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where Nt is the total (or average) number of particles 
employed in the interferometric process. Yet, this is not 
the ultimate limit imposed by Quantum Mechanics. 

In the last few years, it has become clear that quantum 
entanglement has the potential to revolutionize interfer- 
ometry by allowing phase estimations with sensitivities 
up to the Heisenberg limit 0, [1] 
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Pioneering works along this direction were initiated in 
the early 80s 0, inaugurating a large body of litera- 
ture proposing new quantum states and strategies for sub 
shot-noise performances d, @ . Quite recently, several ef- 
forts have been directed to the experimental realization 
of NOON states @, 1, [1 (often called Schrodinger cats 
PUSHES] m the context of trapped ions): 
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The state \ N, 0) contains N particles in mode a and par- 
ticles in mode b (vice versa the state |0, N)), while <p is an 



arbitrary phase. It is widely believed that interferometry 
with the state Eq. ([3]) can estimate unknown phase shifts 
with sensitivity at the Heisenberg limit Eq.@[13]. This 
claim is often accompanied by a simple example. The 
phase shift, induced by an external classical perturbing 
field, is created by the unitary operator Ug = e~~ lJ * e , 
where the generator of the unknown phase translation 
6 is the two-mode relative number of particles opera- 
tor, J z = (N a — Nb)/2. The projection of the new state 
\V N ($)) = U 9 \^n) = (|iV,0) + e* 0If +*)\Q, N))/V2 over 
the initial one gives 



\(^n\^n(0))\ 



2 {N0/2). 
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Orthogonality, (^n\^n(9)) = 0, is first reached at 9 — 
ztw/N, which would suggest that the smallest measur- 
able phase shift is of the order of l/N as well. There 
is a problem, though: in interferometry the incremental 
phase shift, albeit supposedly small, is unknown and the 
phase estimation based on the projective measurement 
Eq.® is ambiguous. Indeed, ($ N \V N (9)) = when 
6 = ±(2n+l)7r/iV, with n = 0, 1, 2, . . . ,N- 1. Orthogo- 
nality alone is not sufficient to determine n, with unpleas- 
ant consequences when trying to estimate the unknown 
value of 9 with an arbitrary large number of particles and 
complete prior ignorance. The 2ir/N oscillation period of 
Eq.((4]) is typical in quantum enhanced technology with 
state Eq.Q. For instance, the multipeak structure of 
Eq.Q is present (even if not generally discussed) when 
measuring the mean value of the parity operator of one 
of the output modes obtained after the state ([3]) has been 
shifted in phase and passed through a 50/50 beam splitter 
0, In this Letter we propose i) a measurement strat- 
egy for the unambiguous estimation of phase shifts with 
uncertainty ~ 1/Nt by using the state Eq.Q, within ii) 
a rigorous Bayesian analysis of the measurement results 
which can be implemented experimentally incorporating 
decoherence and classical noise and iii) maximum priori 
ignorance about the phase shift: — n < 9 < 7r (but we will 
also consider the case of an arbitrary smaller prior). The 
protocol requires p independent interferometric measure- 
ments performed with NOON states having a different 
number of particles, N = 1,2,4,..., 2 P_1 . The sensitivity 
is calculated as a function of the total number of particles 
used in the process, Nt = 2 P — 1. 
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FIG. 1: Phase probability distribution Eq.© with p = 1, 
iVr = 15 (A), and Eq.@ with p = 4, iV T = 15 (B). In 
both cases, the total number of particles is the same, but the 
distribution of B gives a phase sensitivity at the Heisenberg 



by applying the operator e~ t6J '. The final state, ob- 
tained after a further application of Un, is projected over 
\N |). Quantum Mechanics provides the probability to 
have \N |) at the output (which will be denoted as a 
"yes" result), given the unknown phase shift 9 and the 
number of particles, N, of the cat state: P(yes\N,9) 

• N\U N e- l6j *U N \N i)\ 2 = cos 2 (^). Notice that 
this probability is identical to Eq.(|j). The probability 
to obtain a "no" result is, obviously, P(no\N,9) = 1 — 
P(yes|iV, 9). A single interferometric experiment consists 
of p independent measurements. We collect a number p y 
of "yes" and p n = p — p y of "no" results with probability 
Pp(p y ,p n \NT, 9), Nt being the total number of particles 
used in the p measurements. According to the Bayes the- 
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knowledge and P(NT,p yi p n ) provides the normalization 
limit. In C we plot the terms cos 2 (2 fc </>/2) ofEq.®. The solid of P p (4>\NT,Py,p n ), which is a phase probability distri- 



we have P p (<j)\N T ,py,p n )P{N T ,p y ,p n ) = 
y,p n \NT,(p)P(<j)), where P(4>) is fixed by the prior 



blue line is for k = 0, the dashed red line for k — 1 and the 
dot-dashed green line for k = 2. By multiplying these three 
distributions, as in Eq.©, all peaks, except the one centered 
around the true value of the phase shift, 8 = 0, are washed 
out to give the phase distribution D (Nt — 7 and p — 3). 
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From the experimental point of view, the demonstra- 
tion of the Heisenberg limit Eq. requires the creation 
of Schrodinger cat states with a minimum of TV = 8 parti- 
cles, which is within the reach of the present state-of-the- 
art. Cat states with up to 6 ions [12| and 5 photons @ 
have been recently reported. As far as realistic techno- 
logical applications are concerned, however, Heisenberg 
limited interferometry with NOON states Eq.Q would 
likely never overcome the performances of classical in- 
terferometry Eq.([T]), where the typical number of parti- 
cles Nt can be several orders of magnitude larger. We 
therefore extend the previous protocol to reach sub shot- 
noise sensitivity AO ssn /AO srl ~ l/\/2 p — 1, which can 
be implemented experimentally with an arbitrarily large 
number of particles. We address the possibility to reach 
over 3 db sub shot-noise in realistic ion and photon ex- 
periments within the current technology [1, [l(| HH, Q3 , 
both in presence of a strong priori constraint and in the 
more general case of a complete prior ignorance. Our 
results do not only apply to ultra-sensitive interferom- 
etry, but naturally extend to quantum position ing 
clock synchronization [l7|, frequency standards [lj| and 
quantum metrology 

Bayesian analysis with Schrodinger Cats. In the 
following, we discuss the Bayesian phase estimation strat- 
egy considering the Schrodinger cat states realized with 
trapped ions by Wineland and collaborators [l(| [HI, G3 • 
The state \V N ) = (\N |> + i^ +N+1 \N T))/V2 was cre- 
ated by applying the "nonlinear beam splitter" operator 
Un = e 1 ^ Jx e l 5 J ^ to spin-down ions \N I) = | J. 
)i— I 1)ni with f = when N is even, and £ = 1 when 
N is odd [3. The \^n) state is then shifted in phase 
by an unknown quantity 9 (which has to be determined) 



where Nt = Yjj=i Nj and Tj = yes (no) if, in the 
j — th measurement, done with Nj particles, we obtain 
a "yes" ("no") result. Equation© contains all the avail- 
able information needed to estimate 9. We can choose, as 
estimator O es t , the maximum of the phase distribution, 
and, as uncertainty AO, the 68%-confidence interval |6|, 
namely the phase interval containing 68.27% probabil- 
ity given by /q"^ d<j>P p (<f>\N T ,Py,Pn) = 0.6827. In 
the following, for analytical simplicity, we will consider, 
unless explicitly specified, the case in which the mea- 
surement results are only "yes": p v — p, p n — 0. That 
happens with certainty at 9 = 0. The extension to an ar- 
bitrary value of 9, where both "yes" and "no" are acces- 
sible, will be discussed in [2(| , also including decoherence 
effects. 

1 /Nt periodicity of the Bayesian phase distribu- 
tion. Let us first analyze the interferometric experiment 
with a state I^Afy) of Nt particles and a single measure- 
ment: p = 1. The phase distribution becomes 

Pi(0|A T ,p y = l,p n = 0) a cos 2 (N T cf)/2). (6) 

This probability contains Nt peaks separated by a dis- 
tance 2it/Nt, see Fig.(p]A). Hence, our best guess about 
the real phase shift 9 is O est = j^n ± ^ with n = 

0,1,..., Ay — 1, where the error \/2/Nt is the mean 
square fluctuation around a single peak. In practice, we 
do not have any alternative but to choose, as phase esti- 
mator, one of the Nt equivalent peaks of the distribution 
[cf. discussion after Eq.lgJ)]. In this case, however, the 
interferometric experiment does not give any substan- 
tial improvement in phase sensitivity. Tautologically, the 
phase uncertainty would scale with the total number of 
atoms as ~ 1 /Nt only if we have the a priori knowledge 
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that the phase lies in an interval of width 2tt/Nt around 
the real value. The Ay-peaks structure of Fig.([TJA) does 
not allow the Heisenberg limit, even in the case of arbi- 
trary small incremental phase shifts. How is, therefore, 
possible to to select the "right" peak, so to speak, in 
order to build an unambiguous phase estimator? 

The AO ~ N- 3 / 4 limit. In this section we consider 
multiple independent measurements with different num- 
ber of particles. The first measurement is done with 
a cat state of a single particle, N = 1 (with a prior 
knowledge of the phase shift in the interval [— tt,tt]). We 
then perform a second measurement with N = 2, and 
we multiply the resulting distribution with the previous 
one. We repeat this procedure p times, increasing, in 
each shot, the number of particles in an arithmetic se- 
quence N = 1,2, 3, p. The total number of particles is 
Nt = p{p + l)/2, and the phase distribution is 



P p ((/)\N T ,p y =p,p n = 0) OC Y[ cos2 ( -7T 
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(7) 



In the limit of large p, a Gaussian approximation 
gives nLiCOS 2 (M « exp[- f (p+1 f +1) ], where 
YJk=i k2 = P(P + l )( 2 P + i)/ 6 - For p > 1 we have 
Nt ~ p 2 /2, and AO = (9/2) 1 / 4 /N T /4 , which is in good 
agreement with the numerical calculation: 
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This argument can be generalized to the case of any prior 
phase knowledge [— ir/L, tt/L], with L > 1. The goal 
is to obtain an unbiased estimate of 9 within a region 
2n/ L, the peaks outside this region being wiped out but 
the prior knowledge. As before, the protocol requires p 
measurements: the first one with N k, L particles, the 
second one with 2N , the last one with pN . The to- 
tal number of particles is Np(p + l)/2. In the limit of 
large p, we obtain an unbiased phase estimate with a 
sensitivity A6 = With an arbitrary value of 

the true phase shift 9^0, the same protocol provides 
a distribution peaked about 9 with a sensitivity scaling 
as ~ 1/N^/ 4 for Nt 1. This because the cos and 
sin distributions corresponding to "yes" and "no" results 
overlap out of phase and cancel out each other except 
in a region around 9. Yet, even if the scaling ~ 1/N T ' 4 
overcomes the shot-noise Eq.([T]), still this is not the fun- 
damental limit. 

The Heisenberg A9 ~ N' 1 limit. Let us now con- 
sider p independent measurements done with a geometric 
sequence of particles, A^ = 1, 2, 4, 8, ...2 P_1 . The phase 
distribution becomes 

p- 1 /2 k <h\ 

P p ^\N T ,Py=P,Pn = 0) (X J] COS 2 {-f), (9) 

k=0 

where the total number of particles employed in the com- 
plete process is Nt = 2 P — 1. Fig.fTjB) shows the case 



P = 4, Nt = 15, with a prior — ir < 9 < ir. Again, the 
width of the distribution can be simply calculated with 

a Gaussian approximation of each cos 2 ( 2 ^) term, giv- 

ing P p (f\N T ,p y =P,Pn = 0) c expt- ^+y- 1 ) ]. In 
the limit of large p we obtain a phase sensitivity at the 
Heisenberg limit AO = vE/Nt- The numerical calcula- 
tion gives, asymptotically in the number of measurements 
P, 
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for a 68% confidence [21j. We therefore conclude that it 
is possible to obtain an unbiased phase estimation at the 
Heisenberg limit, with repeated measurements and com- 
plete prior ignorance. The trick is to carefully choose the 
number of particles in each measurement and to multi- 
ply the corresponding phase probabilities so to cancel out 
the extra peaks of the phase distribution, see Figs.JI]). 
To clarify this effect, let us consider a phase distribution 
obtained with 2 k particles: cos 2 (2 fc (/)/2). This has max- 
ima in (j) k max {n) = 2irn/2 k , with n = 0, ±1, ±(2 fe - 1). 
Conversely, the phase distribution obtained with 2 fe_1 
particles has minima in 4'^An( n ) = 27m/2 fe with n — 
±1, ±3, ±(2 fe_1 — 1). By multiplying the two distribu- 
tions, the maxima 0^ aa ,(n) superimpose with the min- 
= ±1,±3, .... When taking into 
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account also the distributions with 2 



fc-2 



parti- 



cles, we obtain the cancellation of all peaks 4>max{ n ) with 
n = ±1, ±(2 fe — 1) except the central one, (t>max( n = 0), 
which is enhanced and has a width scaling as ~ l/2 fc ~ 
1/Nt, see Figs.([Tl C,D). This shows that the protocol em- 
ploying A^ = 1,2,4,...,2 P particles is optimal and gives 
the Heisenberg limit with the best prefactor. Indeed, let 



us consider a general sequence A^ = l,r,r 



rP~\ with 



an arbitrary integer r. If r = 1, we recover the shot- 
noise limit Eq.JT]), for r — 2, we have, as already dis- 
cussed, a perfect superposition of maxima and minima. 
For r > 2, the phase distribution is characterized by sev- 
eral, strongly weighted, peaks which increase the phase 
uncertainty. To overcome this problem we must repeat 
AI > 1 times each interferometric measurement with a 
NOON state of fixed A^ before increasing the number of 
particles (employing a total Nt — M X)fe=o rk ) • F° r a 
sufficient large M, it is possible to recover the Heisen- 
berg limit Eq.([2]), but at the price of a higher prefactor. 
In fact, we have that A9 ~ vM/Nt, as a consequence 
of the statistics of independent measurements. In Fig.@ 
we analyzed the sensitivity obtained with different r as a 
function of M: for r = 2 (black circles) the best perfor- 
mance is obtained at M = 1, for r — 3 (blue points) at 
M = 4, for r — 4 (red diamonds) at M — 6 and for r = 5 
(green squares) at M = 9. 

As mentioned before, the analysis can be extended to 
estimate an arbitrary unknown phase shift 9^0. How- 
ever, in contrast to the AO ~ N^ 3 ^ 4 case, here we have 
to consider multiple repeated measurements in order to 
concentrate the probability in a single peak, even for the 
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FIG. 2: Phase sensitivity obtained with iV = l,r, r 2 , r- p_1 , 
where each measurement (using N = r k particles) is repeated 
M times. The best strategy is obtained when r — 2 (black 
circles) and M = 1, due to a perfect superposition of max- 
ima and minima in the Bayesian distributions, see Fig.(JT]D). 
Notice that the optimal value of M, as well as the phase uncer- 
tainty, increases with r: here we consider r = 3 (blue points), 
r — 4 (red diamonds) and r = 5 (green squares). Lines are 
guide to the eye. 



case r — 2. While the scaling ~ 1/Nt is preserved [20| . 
Eq. (JTUJ) gives a lower bound of phase sensitivity. 

Sub shot-noise with Beryllium ions. The experi- 
mental demonstration of the Heisenberg limit would re- 
quire the creation of Schrodinger cats having N = 8 par- 
ticles. The biggest Schrodinger cat created experimen- 
tally so far is with N = 6 9 Be + ions. This is still suffi- 
ciently large to reach a sub shot-noise phase sensitivity 
[22| . In the following, we demonstrate, using the fidelities 
measured experimentally in [Tol . [TP . [l2j with Beryllium 
ions, the possibility to reach a phase sensitivity gain of 
0.8 db with respect to shot-noise for a priori — 7r < 9 < tt. 
The protocol is quite similar to the one discussed above 
and requires Bayesian probabilities calculated with dif- 
ferent number of ions combined with M replica of the 
measurement process. We remark, however, that, in 
the analysis, we should now replace the ideal probabil- 
ity distributions with those retrieved experimentally. In- 
deed, we need to include the effects of noise and deco- 
herence present in the experiments which will inevitably 
decrease the sensitivity of the interferometer. This step 
can be considered as a "calibration" : once the experimen- 
tal probabilities are retrieved and inverted with Bayes, 
the interferometer is ready for its use. Many sources 
of experimental imperfections and decoherence conspire 
against phase uncertainty enhancement with cat states. 
Laser intensity fluctuations and magnetic field noise have 
been discussed in [T^]. Imperfections in the creation of 
the state Q due to non ideal U n and decoherence have 
the common effect to decrease the oscillation amplitude 
of P(yes|iV, 9) [24]. A fit of the experimental probabil- 
ities with P exp (ycs\N, 9) = A + C N ^ Nl /2cos(N9) and 
P ex Jno\N 9) = 1 - P exp (yes\N,9) has been given in 
fTfjl. [TTl [T2T ] with the values of A and Cnt,ni reported 
in [25|]. The main effect of the experimental noise is to 
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FIG. 3: Phase probability distributions obtained by combin- 
ing M times (in (A) M = 1, in (B) M = 10) the Bayesian 
distributions for N = 1, 2, 3, 4, 5, 6 particles and 9 = 7r/2. 
The dashed red line is the ideal case, while the solid blue line 
has been obtained with the experimental fidelities, see (25|. 



decrease the fringes contrast Cjv|,jv|j which, in the ideal 
case, is equal to one. The experimental Bayesian phase 
probability distribution associated to a detection of a 
"yes" or "no" result, P exp (cf)\N, {yes, no}), is obtained in- 
verting P({yes, no}| N, 9). We are now ready to simulate 
a realistic phase estimation experiment with Beryllium 
ions: i) A "yes" or "no" result is chosen with probability 
P({yes, no}|7V, 9) with an unknown, but fixed, value of 9 
and for different number of ions, N = 1, 2, 3, 4, 5, 6; 
ii) We repeat these measurements M times; iii) We cal- 
culate the Bayesian distribution associated with each 
"yes" /"no" result and values of N; iv) We multiply all 
Bayesian distributions obtained in iii). This provides the 
final phase probability, from which we retrieve the esti- 
mated value of the phase shift and its confidence. The to- 
tal number of particles used in this process is Nt = N p AI, 
with N p = Yl=i k = 21 - Ideally, with N = 1, 2, ...,p, the 
phase sensitivity would scale as 



A6 = 



(9/8) 



1/4 



(9/8) 
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(11) 



where N p = p(p + l)/2. In Fig. ([3]) we plot the theo- 
retical and experimental phase probabilities for M = 1 
and M = 10. Notice that the experimental distribution 
for M = 1 is characterized by a large tail. This arises 
from the reduced fringe visibility and it would strongly 
dilute the confidence of the phase estimation. On the 
other hand, by multiplying the M — 10 distributions, 
step ii), we strongly decrease the weight of the tail with 
respect to the central peak, see Fig.(j3lB), at the price, of 
course, to increase the shot-noise. It is worth to empha- 
size, however, that, because of noise and decoherence, at 
the end of the day we could, in principle, get an exper- 
imental phase sensitivity even worse than the classical 
shot noise. Asymptotically in M, the phase probability 
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Eq.© can be written as 



Pex P (ct>\N T ,0) = n Pex P (cj>\N,yes) 



MP emp (y CS \N.e) 



xP exp (^no) MP -( nc W). (12 ) 

In this limit, and with ideal fidelities, we would have a 
phase independent gain Gth — 3.18 db with respect to 
the shot-noise limit AQ sn . Such a gain would be compa- 
rable with the best performances obtained to date with 
photons [2(1 . As a consequence of imperfections and de- 
coherence, the experimental gain is lower than the ideal 
prediction and depends on the phase shift, see the solid 

0.75 



exp 



red line in Fig. - The maximum gain is G. 
dB, around the optimal working point 9 ~ 0.3. In princi- 
ple, an even higher gain can be obtained using states with 
N = 1,2,4,8, ...,2 P ~ 1 particles: the sensitivity would be 
bounded by the ideal value 



AO = 



1. 
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with N P = 2P- 

So far we have considered the case of complete prior 
ignorance. While this can be important for technological 
and basic science applications like in gyroscopes or clock 
synchronizations, it is sometimes possible to confine the 
priori within an interval —ir/L < 9 < tt/L, with L > 1. 
In this case an unbiased phase estimation can be obtained 
with M replica of the interferometric measurement, each 
with a Schrodinger cat state of a fixed number of parti- 
cles N ~ L. With ideal distributions, we would obtain 
a phase-independent sensitivity AO = with 

Nt = M N. The gain with respect to the shot-noise, 
obtained with the experimental probabilities, is shown in 
Fig. (U)) where the arrows indicate the prior knowledge 
tt/N for the various cases. The gain is maximum at 
9 = 7r/ (27V) where the experimental probabilities to have 
a "yes" or "no" result are closer to the ideal ones. No- 
tice that the sensitivity at first increases and eventually 
decreases with the number of ions. This is caused by the 
competition between the gain obtained increasing N and 
the loss of visibility due to decoherence, which is higher 
for larger cats. With the experimental fidelities measured 
so far, the best scenario is obtained employing cat states 
of N — 3 particles, which provides a gain up to 1.63 db 
with respect to the shot-noise. Yet, it is clear that, with 
a priori —n/L < 9 < n/L and with cat states having a 
number of particles N ~ 2 k L, with k = 1, 2, 3, we can, 
in principle, further decrease the phase uncertainty than 
using cats with a fixed number of particles. 

We remark that the total number of particles used in 
our phase estimation protocol can be arbitrarily large 
(being proportional to the number of replica M). This 
is important for realistic technological applications since 
the number of particles in Schrodinger cat states would 
probably remain relatively small, at least in the next fu- 
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FIG. 4: Gain (db) with respect to the shot-noise limit 
obtained with the experimental parameters [25( ] . The blue 
dashed line is the case N = 2, the green dot-dashed line 
N = 3, the brown dot-dot-dashed line N = 4, the sky-blue 
dot-dot-dot-dashed line N — 5, and the black dotted line 
N = 6. Colored arrows indicate the upper bound to the phase 
prior for each data set. With maximum priori ignorance, we 
combine M replica of the Bayesian phase distributions with 
N = 1,2,3,4,5,6, Eq.JTJJ, (solid red line). The horizontal 
dotted line gives the shot-noise. 



ture. We can expect to have, in a few years, the produc- 
tion of robust high-fidelity cat states up to ~ 10 photons 
or ions, allowing the saturation, at least in principle, of 
Eq.([2]). Major obstacles to these efforts are creation im- 
perfections and decoherence 0, H3] • On the other hand, 
Bose Einstein Condensation might offer the possibility 
to create NOON states with a larger number of parti- 
cles [1^]. A different, promising strategy to experimen- 
tally reach sub shot-noise sensitivities is to use number- 
squeezed states, which have been recently experimentally 
demonstrated [29l . [30j with up to few hundred neutral 
atoms. In any case, crucial to the success Heisenberg 
limited interferometry, is the realization of highly effi- 
cient number counting detectors, which can be prob ably 
developed in the next generation of experiments [30ll3ll|. 

Conclusions. Ultra-sensitive interferometry requires 
unbiased phase estimation protocols and carefully en- 
gineered maximal quantum correlations among input 
states. In this paper we have developed a novel 
Bayesian protocol based on multiple measurements with 
Schrodinger cat states of variable number of particles. 
This achieves the Heisenberg sensitivity AO ~ N^ 1 with 
an arbitrary prior phase uncertainty. This protocol over- 
comes difficulties present in previous approaches where 
the estimate was strongly ambiguous, so to require a 
prior knowledge of the true value of the phase in the 
restricted interval —tt/Nt < < tt/Nt- We have also 
demonstrated the possibility to obtain sub shot-noise 
phase sensitivity with the experimental data published in 
[ToL [Til OjJ on the creation of Beryllium ions Schrodinger 
cat states. Our results do not depend on how the state 
Eq.Q is created, nor on the specific interferometric ap- 
paratus as far as the phase probability distributions have 
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an oscillating pattern with a period depending on the 
number of particles. 
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